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A wide range of equations related to free surface motion in two dimensions exhibit the formation 
of cusp singularities either in time, or as function of a parameter. We review a number of specific 
examples, relating in particular to fluid flow and to wave motion, and show that they exhibit one 
of two types of singularity: cusp or swallowtail. This results in a universal scaling form of the 
singularity, and permits a tentative classification. 
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^ ' I. INTRODUCTION 
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o , 

The geometry of singularities of wavefronts, i.e. caustics, has been much studied nj. The caustics, themselves, 
exhibit singularities such as cusps, which can be classified using catastrophe theory [2j, l3j. In the present review we 
point out that a similar classification appears to apply in a variety of nonlinear PDE problems, which describe the 
motion of a free (fluid) surface. 

The two-dimensional problems studied here can be written as a mapping of the physical domain to the unit disc, the 
. free surface being represented as the circle. The appearance of a singularity is associated with the non-invertibility of 
this conformal map at a given value of time or of some other parameter. Since the problems studied here are endowed 
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with a holomorphic structure, such non-invertibility has a very particular geometrical consequence: the formation of 
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a cusp with a 2/3 exponent or of a swallowtail singularity. When looking at the mapping as it evolves in time, or 
depends on a parameter, it is apparent that the mapping remains smooth as the self-intersection occurs. For this 
reason the curve can be written as an analytic function of a suitable parameter. Performing a local expansion around 
the point of self-intersection, the properties of the curve are determined by the lowest order terms of the Taylor 
7-H ■ expansion in the parameter. 

^ ' The simplest case of self-intersection is that of the cusp. Let us assume that the coordinates x and y can be 
, expanded into a Taylor series as function of some parameter, 6, and that we are interested in the neighborhood of 
9 = 0. Since we describe phenomena up to arbitrary translations, the generic description is x = a\0 and y = a^Q. 
However, by performing a rotation one can always make sure that one of the coefficients (a2, say), is zero. The 
singular case corresponds to the situation where a\ vanishes as well, so we put x = eO, with the singularity at e = 0. 

Of course we are now required to expand to higher order. The next non-trivial order in y gives y = 2 /2, where 
the coefficient can be normalized by rescaling 9. The expansion in x has to be pursued to third order, otherwise the 
curve would be degenerate for e = 0. Thus we finally have: 
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x = e6 + a6 3 /3, y = 6 2 /2, (1) 
where a is an arbitrary constant. The quadratic coefficient in the expansion of x can be eliminated by a shift in 



(a) (b) (c) 

FIG. 1: The formation of a cusp, as described by {l}. Shown is a non-intersecting curve (e = 2), a 2/3 cusp (e = 0), and a loop 
(e = —1), from left to right. 
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6, with subsequent rotation and translation. Thus up to translations and rotations, {T]) is the only generic way the 
self-intersection of a curve in the plane can occur, as illustrated in Fig. [TJ For e < the curve self-intersects, at 
the critical point e = a cusp is formed. It is clear from ([I]) that this cusp has the form y = (x/a) 2 ^ 3 , i.e. it is 
associated with a universal power law exponent. In catastrophe theory ||, {!]) is the curve associated with the "cusp" 
catastrophe. 

As we will see below, higher order singularities can also occur, which presumably is due to a hidden symmetry of 
the problem. As a result, the coefficient in front of the quadratic term vanishes in (JTJ) : y — e s 6 2 /2 + bx6 3 /3 + b 2 9 A /4. 
Once more, the cubic term can be made to vanish by performing a rotation, so we arrive at the following generic form: 

x = e s 9 + a0 3 /3, y = e s 9 2 /2 + b9 4 /4, (2) 

which exhibits a much milder singularity. Namely, for e s = the behavior is like y oc x 4 / 3 . For e s < 0, the curve splits 
into two cusp singularities of the form shown in Fig.[TJ where the case b > a corresponds to e > in the representation 
([T]), b < a to e < 0. In all the cases to be discussed below, we always find the particular case a = b, which means for 
e s < the solution is exactly at the cusp singularity. This is the swallowtail known from catastrophe theory, which 
can be seen as a collision of two cusps (see Fig. [2]). It has the universal form 

x = e6 + a9 3 /3, y = e9 2 /2 + a6» 4 /4, (3) 

once more with e, a being parameters. 

The two catastrophes (jTJ) and ([3]) can occur either as function of time or of some control parameter. Let us introduce 
the time distance t' = to — t to the singularity, and assume the scaling |e| = |i'| 7 . The time before the singularity 
t' > corresponds to e > 0, and vice versa. Then 

x=\t>\^/ 2 X Cl y=\f\"<Y C) (4) 
X c = ±<r + ao- 3 /3, Y c = (7 2 /2, (5) 

is the self-similar form of the cusp singularity, and 

x=\t'\ 3 ^ 2 X s , y=\t'\ 2 ->Y s , (6) 
X s = ±g + a<7 3 /3, Y s = ±<r 2 /2 + acr 4 /4 (7) 

of the swallowtail singularity. In each case, the + or - signs correspond to the similarity function before or after the 
singularity, respectively. We remark that the cusp singularity can be written as the following cubic: 

X 2 = 2Y c {l±2aY c /3) 2 . (8) 

As we have mentioned before, the similarity function ([B]) valid after the singularity contains two cusp singularities, 
which occur for a = ±1/ ^fa. Thus if one writes a = ±1/ ^Ja + s, shifts the cusp to the origin and performs a rotation, 
one obtains to leading order in s: 



/ ±as 3 \_f 1 T\A \ f X s ± 2/(375) \ 
\(a + l)s 2 )-{±V^ 1 J\Y S + l/(4o) J 



(9) 
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Equation 


Type 


7 


Section 


Wave fronts 


swallowtail 


1 


m 


Hele-Shaw flow 


cusp 


1/2 


El 


Potential flow with free surface 


swallowtail 


// 


w 


Porous medium equation 


cusp 


// 




Viscous flow with free surface 


cusp 


// 


m 


Born Infeld equation 


swallowtail 


1 


vm 



TABLE I: A summary of evolution equations discussed in this paper. The classification as "cusp" or "swallowtail" refers to |T|) 
or @, respectively. In the case of time-dependent problems, the exponent 7 is defined by Q or (|6}, depending on the type of 
singularity. 



The + or - signs correspond the left and right cusp, respectively. This demonstrates that (X S ,Y S ) locally traces out 
a cusp after the swallowtail singularity, as seen in Fig. O 

Table [J summarizes the the problems to be studied in this paper, and cites the relevant sections. Each equation, to 
be discussed in more detail in the sections below, exhibits singularities which can be classified as either being of the 
"cusp" or the "swallowtail" type. Some are evolution equations, in which case we give the temporal scaling exponent 
7, others exhibit singularities as function of a parameter. In each case, we attempt to give physically motivated 
examples. 



II. THE EIKONAL EQUATION 



A. Hamilton-Jacobi equation 



We consider the simplest case of wave propagation in a homogeneous medium, so the eikonal equation becomes 
V n = 1: the normal velocity of the wave front is constant everywhere. We will only consider wave propagation in two 
dimensions, so the wave front is a curve. A straightforward way to solve the eikonal equation is to consider y — h(x, t) 
so that one can write the equation in the form 



ht = VTThJ, h(x, 0) = h (x), (10) 

which expresses the condition of constant normal velocity. 

To solve (fTUj) . we pass to the "particle" description; the equation for the wave front (fTOjl is the Hamilton-Jacobi 
equation for the Hamiltonian 

H(x,p)=-y/l+p>, (11) 
where p — h x . The initial condition (xq, h' Q (xo)) yields a curve (x,p) in phase space, and the trajectories will be 

dx dH p dp dH 



dt dp y/i 4- p 2 ' dt dx 
Thus the solution to the particle problem is 



= 0. (12) 



p = const = h'n(xn), X = Xn ; = t — Xn o( o) ^ 

° l 7 y/TT]? y/1 + WW) 



We can obtain an explicit solution by noting that 



dh dh dxo ^, . > 
dx dxo dx 



where from (|13[) 

dx y/l + fr3 



dx = i _ fep f 



FIG. 3: A caustic in a coffee cup. On the left, physical reality. On the right, the evolution of the wave front, as it passes 
through the swallowtail singularity. The caustic forms a 2/3-cusp, which corresponds to a swallowtail in terms of the shape of 
the wavefront. 



Integrating the resulting expression for dh/dxo, we finally obtain 

X=X °--^W^ h=hM+ vm§- (14) 

This is an explicit solution of (flU)) , which has the additional advantage that it can be continued across any singularity 
the wave front may encounter. 

A caustic is a place where rays meet, and thus dx/dxo = 0, dy/dxo — 0. The two conditions turn out to be 
equivalent, and one obtains 



to = - (15) 

"0 

It is also clear that points on the caustic correspond to singularities of the wave front [ij . The first singularity occurs 
at a time t corresponding to the maximum of the curvature, for which there is optimal focusing. The universal 
structure is obtained by expanding ho in a series, which without loss of generality only contains even terms: 

ho = aix'l + ci 2 Xq + (16) 

The first singularity occurs for to = l/(2ai) + . . ., and an expansion leads precisely to ((5]), namely 

x = t'a + aa 3 /3 (17) 
y = l/(2ai) -t' + t'a 2 /2 + aa 4 /A, 

where a = 3(af — a2)/(4af), and a — 2a\Xo] the scaling exponent is 7 = 1. After the singularity, t' < 0, (fT7|) is a 
swallowtail, and the two cusp points trace out the caustic. This corresponds to the positions o~ c — ±y-t7a, and so 
to leading order the equation of the caustic is 



_ 4i /3/2 x 
± — > ^=2^-*' (18) 



which is a normal 2/3-cusp, see Fig. [31 



B. conformal mapping 



There is a second way of solving the eikonal equation, closer to the methods to be used in the next sections, based 
in the use of a complex representation of the front. We present it here to stress the similarities between the different 
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problems studied in this paper. If one identifies the points of the curve (x(a, t),y(a, t)) as points z(a,t) — x(a,t)+ 
iy{a, t) in the complex plane, and the velocity u = (u, v) as the complex number u + iv, then 

z t — u + iv. (19) 

The tangent and normal vectors have a complex representation 

t = z s , n = iz s , (20) 

where s is the arclength parameter. Hence we can write 

u • n = Re (iz s zi) = -Im (z s zt) , 

and the eikonal equation is simply 

lm(z s zi) = -l, (21) 

with the constraint |z s | = 1. Equivalently, one can introduce an arbitrary parameter a of the curve, and (|2ip becomes 
or, equivalently, 

ds 

lm(z a - t ) = - — , (22) 
aa 

without the necessity of introducing a constraint. 

We now demonstrate that ([22]) yields the same local solution as before, by verifying that (fTTf yields, at leading 
order, a solution of (f2"2"|) : 

z a = (t' + aa 2 ) + i(t'a + acr 3 ) , z t = -a + i(l - cr 2 /2) 

2 

Im(z CT zT) = -(1 + y)(f' + aa 2 ) 
ds 



— = y/(t' + OfT 2 ) 2 + <r 2 (t' + aa 2 ) 2 = (1 + —)(t' + aa 2 + 0{a 4 )). 
da 2 

III. HELE-SHAW FLOW 

A Hele-Shaw cell consists of two closely spaced glass plates, partially filled with a viscous fluid. The problem is to 
find the time evolution of the free interface between fluid and gas. Here we consider the case that the fluid occupies 
a closed two-dimensional domain Q. Within Q, the pressure obeys Ap — 0, with boundary conditions 

p = (23) 
V n = -Vp-n (24) 

on the free surface <9f2. We write z = x + yi and 

p = Re$(z) 

together with the conformal mapping 

z = f(.tt), t = re w , (25) 

which maps |£| = 1 onto dfl. If we consider 

$(z(£, t)) = log^ inside O, 

then condition (|23[) is automatically satisfied. Moreover 

fd® \ (d§ 1 \ / 1 . \ 

— \p ■ n = —Re — — iz s = —Re —— — iz s = —Re - — iz s , 



dz J \ d£_ z^ / \^z^ 

where s is the arclength parameter. Notice that one can write for |^| = 1, using £ = e 18 ^'^ 
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and hence 

- Vp • n =e B . (26) 

Since 

V„ = Re (ziizs) = -O s Re 
we arrive, combining this with (|26p . at the equation 

Re(z^) = -1. (27) 

Noting that £z£ = iz#, we finally obtain 

Im(z e zT) = -l. (28) 



Equation (|28|) is identical to (|2Tj) except for the fact that the space variable is not the arclength parameter s, but the 
an arbitrary parameter 9. The non-invertibility of /(£, f) signals the appearance of a singularity. 

Now we will present an exact solution of (|2T|) . using a particular (polynomial) ansatz for the mapping [231 H 
We consider the simplest case of a quadratic: 

f(^t)^ ai (t)^ + a 2 (t)e, (29) 

and show that it leads to cusped solutions. However, we expect cusp formation to be a generic feature. The reason 
is that the formation of a singularity is associated to the non-invertibility of the conformal map /(£, t) and this is 
equivalent to /'(£, t) being zero at some point £ (at the time t of formation of the singularity), which leads to a 
generic quadratic behavior of /(£, to) = &o + b 2 (£ — Co) 2 near £o- Other local expansions of /(£, to) around £o ; of the 
form /(£,to) = b + fr„(£ — £ ) n with n > 2 are also possible. They would lead to different form of cusp, but cannot 
be generic since small perturbations of the initial data would produce quadratic terms. 
Inserting ([23)) . (f2"9"|) into (f2"8"|) , we find a solution if the following system of ODEs is verified: 



aia\ + 2a 2 a' 2 = -1 (30) 
axa 2 + 2a 2 a[ = 0. (31) 



Direct integration of the equations leads to 



2 



2 

A singularity occurs when / fails to be invertible, that is 



l a\+ I L=A-t (32) 
A = \a\{0) + o2(0) , B = a 2 (0)a 2 (0) (33) 



or equivalently when 



f = a 1 + 2a a £ = 0=K = -^- = -l (34) 
at, 2a 2 



ai -4^ =0^ ai = (25)5,02 - i(2B)s. (35) 



From (|32[) one can compute the singularity time 

to=A- |(2B)i = ia 2 (0) + a 2 (0) - ^(2a 2 (0) a 2 (0))i. 
Now let us consider a particular solution, for instance 



so that 



oi(0) = 1 , oa(0) = i (36) 



129 „ I /3\ 4 , , 1 , x 1 



A= 256' S= 16' <0= UJ ' ai(< ° )= 2 ' a ^ = I 
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FIG. 4: The formation of a cusp in Hele-Shaw problem with suction at the origin. 



This implies the formation of a singularity at x(w, to) — — \ and y(7r,<o) = 0. 

Local analysis leads to ai(t) — 1/2 + a(t), where 3a 2 (t) ~ t n — t. Thus the scale factor a{t) is 

o(t)~-L(t -t)i = -L t 'i. (37) 

and a 2 {t) = l/(16a? (t)) ~ 1/4 - a(i). Writing £ = -1 + £, we deduce 

/(£, t) = ai(t)£ + a 2 (i)£ 2 ~ -1 - 2o(t) + 3a(i)£ + i?. (38) 



Since 



one has 



-1 + e = -e 1 ^ = -1 + (1 - e 49 ) ~ -1 + y - ^ - i0 + i 9 l + O(0 5 ) 



q2 M _ fl 3 



/(£, i) ~ -1 - 22(t) + 3o(t) [ y - ^ - + I - ^ 2 + ^ - \0\ + 0(6 5 ), 

leading to the leading order contributions (all other terms are subdominant for small values of a(t) and 6): 

x(9,t) + 1 25(f) - -9 2 (39) 

y(9,t) ~ -35(f)0-i# 3 , (40) 

which is the desired local solution we have been looking for. 
Defining 

Y _ y_ Y _ x + 1 1 o_ f41 s 

one obtains the cusp ((H), with a — 9V3 and the scaling exponent 7 = 1/2. A different choice of initial conditions |36|) 
will of course lead to a different value of the parameter a. In Figure[3]we represent the interface profiles corresponding 
to the example developed above at the initial time and at the time of formation of the cusp. 
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IV. POTENTIAL FLOW 



We consider the two-dimensional flow of an ideal fluid below a free surface. Inside the fluid, the fluid velocity 
u = (u,v) satisfies 

u = V<f>, = 0, (42) 

where <f> is the velocity potential. The free surface is convected by the fluid velocity, and the free surface is at constant 
pressure. We consider the simplest case of steady flow, as well as no body or surface tension forces. According to 
Bernoulli's equation, the fluid speed then has to be constant on the free surface. 

Exact solutions to the flow problem can be found if the fluid domain is bounded by free surfaces and straight solid 
boundaries alone @, by mapping the fluid domain onto the upper half of the complex plane, which we denote by £. 
Here, we consider only the even simpler case of only free boundaries, i.e. that of a two-dimensional fluid drop. To 
this end, one introduces the complex potential w = (fr + itp, where ip is the stream function. The derivative of w gives 
the fluid velocity: 

dw ■ - -ie fA ?\ 

-—=u — iv = qe , (43) 

dz 

where q is the particle speed. 

On the fluid boundary, ip is constant, and this constant can be chosen to vanish. Thus w(C) is real on the real axis 
(the boundary of the fluid drop), and so dw / dC, must be real as well. 

Following Hopkinson [To| we drive a fluid motion inside the drop by placing singularities inside the drop. We will 
consider the case of a vortex dipole and a vortex at the same point. We have 

dw dw dz 

where dz/dC, contains no singularities in the upper Q plane, since the representation of the fluid flow must be conformal. 
But this means that dw/dC, must have the same singularities as dw/dz, except that we now have the freedom to choose 
the position of the two singularities and the orientation of the vortex doublet arbitrarily. Thus if we choose the two 
singularities at £ = i, and the orientation of the doublet toward the positive real axis, the potential for the doublet 
must locally look like w oc — i) — imln(£ — i). Here m is the relative strength of the vortex. 

To insure that dw / dQ also obeys the boundary condition, one has to add image singularities at £ = — i: 



dw I I / 1 I 



dC (C-*) 2 (C + *) 2 "AC-* C + i 

Depending on the value of m, two different cases arise. For simplicity, we only consider the case m < 1, in which 

dw _ 2(1 + m)(C 2 - 7 2 ) 



(45) 



(46) 



Where 7 = — m)/(l + m) real. 

Manifestly, dw/dC, is real on the real axis, and has the right singularities at £ = i. However, the information 
contained in (|45|) is not enough to reconstruct the mapping z(Q we are after. Following Kirchhoff [llj and Planck 
[T^ |. we define another function 

d,z 

il = In— ee -lnq + ie, (47) 
dw 

which will also be represented in the £-plane. Since q is constant along the free surface, and choosing units such that 
5=1, the function Q will be purely imaginary on the real £ axis. To find f2, we once more proceed in two steps. 
First, we find the singularities of Q. From the definition, 

„ , dw , dz , . . 

^ = -m- + m-, (48) 

where the second contribution is conformal in the upper half plane. This means f2 has singularities only for £ = i, 
where it behaves like Q oc 21n(u — i). Second, we have to make sure that is imaginary for real C, which is achieved 
by 

= 21n^^-. (49) 
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FIG. 5: On the left, we show the entire drop for m — 1/3. On the right, a closeup of upper part of the drop for m = 0.94. 
The swallowtail has opened, and there is no more self-intersection. 



Now we use the fact that 

dz dz dw n dw 2(1 + m)(£ 2 — 7 2 ) 

dC"d^dC =e dC = (C + «) 4 ■ { ' 

This expression can be integrated to find the transformation z(£) between the fluid domain and the £-planc: 

z 1 i 1 + 7 2 



2(1 + m) C + i (C + i) 2 3(C + *) 3 ' 
From the real and imaginary part of this expression, we obtain 



(51) 



, = -(2(i +ra ))c5^i!£±ii! (52) 

The mapping is such that C = ioo gets mapped to the origin. A typical drop shape, for m = 1/3, is shown in Fig. [5] 
(left); it exhibits two cusps. This feature is generic, in that it exists for a continuous range of values 0.9427 > m > 0. 
On the right of show a closeup of the top of the drop, close to the upper end of his range. 

The origin of the double cusp lies in a swallowtail transition for m sa 1 that occurs at the point x = Q, y = 4/3 in 
real space. Namely, a local expansion of (|53|) gives with m = 1 + e 

.x = 2eC + 4C 3 /3 (54) 
y = 4eC 2 + 4C 4 , (55) 

where £ is taken as real. The same expression results if the shape of the drop valid for m > 1 is expanded locally. 
The corresponding swallowtail is shown in Fig. [S] using the full solution Fig. 1531 For m very close to the transition, 
the free surface must self-intersect, as shown in Fig. [6] (right). However, below a value of about m w 0.9427 the 
self-intersection disappears. 

Another exact solution of potential flow that exhibits a cusp, but in the presence of gravity was found by Craya and 
Sautreaux [H, Qj]) Gil ; liquid is layered above a two-dimensional ridge. At the crest of the ridge, which has opening 
angle 27r7, there is a sink. For the special case 7=1/3 there is an exact solution, given by 
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FIG. 6: The birth of the swallowtail, as described by (|53[) . For m = 1 (left) one is exactly at the transition, for m = 0.98 
(right) one finds a swallowtail. 



where I = exp(— iff), —it < 9 < it. From a local analysis we find that 

o2/3 o2/3 

x = -—6 3 , y = -^9 2 , (57) 
36 12 ' y ' 

which is once more a 2/3 cusp. Of course this is to be expected, since gravity cannot change the local behavior near 
a cusp. However it is not clear where this comes from in terms of the swallowtail described above. 



V. POROUS MEDIUM EQUATION 

Another problem, closely related to the above, concerns the two-dimensional flow of oil in a porous medium. The 
oil is layered above heavier water, and is withdrawn through a sinkhole. The interface between the oil and the water 
is deformed, and forms a cusp at a critical flow rate. Inside the oil domain, one has to solve Laplace's equation for 
the velocity potential, cf. (|42|) . In the stationary case, one has the usual condition of vanishing normal velocity on 
the free surface. However, from a pressure balance which includes the hydrostatic pressure, one obtains [l6l [l7|: 

u 2 + v 2 = Kv, (58) 

where u, v are the horizontal and vertical components of the velocity, respectively. 

Once more hodograph methods can be applied, but the problem can be solved only in the critical case at which 
there is cusp. Namely, in the subcritical case where there interface is smooth, u = v = at the stagnation point below 
the sink. In the presence of a singularity, u — 0, v = K at the singularity, and the free surface gets mapped onto a 
circular arc in the hodograph plane. The resulting interface shape is 



2 [ " 1 x " 1 + arctanhV^Tl) ( — — ) (59) 



C7T 2 J_g \ 07 / \S- (7; a — a 



y = -(ln^-ln^l±l\, (60 ) 



ctt \ a — a s 1 — j3 — a & 

where a parameterizes the curve and 07 and a s are determined from implicit equations involving c and /3. 
For (3=1 the curve defined by (|5"9" ]) . (|60p develops into a cusp. A local expansion yields 



/ a 2A 
VTT^da = — (cr + 1) 3/2 , 



(61) 
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where A is a constant. The y-coordinate is linear in t + 1, thus (f6"Tj) describes the usual 2/3 cusp. For j3 < 1 the 
interface self-intersects, at the line of symmetry, so this appears to be the generic cusp scenario. Numerical results 
confirm that the curve is smooth before the cusp forms (subcritical case), and a cusp forms in agreement with the 
exact solution ({55)) . (|6T))) . 



VI. VISCOUS FLOW 



Here the flow is governed by the Stokes equation, which in two dimensions can be written in terms of the stream 
function ip, with u = ip y and v — —ip x - The stream function obeys the biharmonic equation A 2 ip — 0. In a stationary 
state, which we are considering, the surface of the fluid is a line with ip = const. On this surface, we also have the 
surface stress condition 

GijUj = JKUi, (62) 

where 7 is surface tension and k the curvature of the interface. Once more the flow is driven by singularities, such as 
a vortex dipole [l^ . 

A complex formulation of this problem was developed by Richardson [191 ] . The stream function is written as 

ijj = Im{f{z) + zg(z)), (63) 
where / and g are analytic functions. The boundary conditions at the free surface can be shown to be 



Im 



f{z) + zg{z) 



1, f(z)+zg(z)=0, (64) 



where r\ is the viscosity of the fluid. 

In [18j , the complex formulation (|64[) was used to calculate the following model problem: a vortex dipole of strength 
a is located at a distance d below a free surface of infinite extend. Surface tension is included in the description, 
but the effect of gravity is neglected. The deformation of the free surface by the viscous flow is determined by the 
capillary number 

Ca = ^L, (65) 

which measures the ratio of viscous forces over surface tension forces. The solution of the problem is too involved to 
be presented here. The exact surface shape, in units of d, is given by the function 

cos 6 

x = a cos 9 + (a + 1) -, (66) 

1 + sin 

y = a(l +sin<9). (67) 
The parameter a is determined from the equation 

l ,Ca= Z^±^m , (68) 
1 + a + ^-2a(a + 1) 



where 



m = 7 , I A , — , 7T77 ( 69 ) 



{-2a/ {a + 1))V4 + (( a + l)/(_2a))V 
and K is the complete elliptic integral of the first kind: 

K{m) = / - (70) 
J a yl — m z $m 6 

In (f6"8")l we have only reported the form of the equation for the more relevant case Ca > 0. Asymptotic analysis of 
([68 ]) -([70 ]) reveals that for large Ca, 

1 32 

a =-- + e, e w — exp{-167rCa} . (71) 
o y 
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FIG. 7: Experimental data on the cusping of a viscous fluid, taken from [2Cj- On the left, a closeup of the tip of a cusp on the 
surface of a viscous fluid; the scale bar corresponds to 200 fi m. On the right, the radius of curvature of the (almost) cusp as 
function of capillary number. In agreement with (|74p . the dependence is exponential. 

It is easy to confirm that (|66[) yields a cusp for a = —1/3, i.e. for Ca = oo or vanishing surface tension. If one 
expands around the cusp point by putting — n/2 + 6, one obtains 



2e 6 3 

x = o 

3 12 

2 S 2 



(72) 
(73) 



which is the generic cusp scenario (p}. As is apparent from Fig. [1] the case e < leads to self-intersection of the free 
surface, which is of course not physical. The radius of curvature of the cusp for e > is given by 



256 

Rfn — exp{-327rCa}, 

as found from (|72 |) . (TfT ]) . The exponential dependence ([74]) has been confirmed experimentally (cf. Fig. [7|). 



(74) 



VII. BORN-INFELD EQUATION 

The ideas presented here are not restricted to free surface problems. As an illustration we present a problem that 
appears in connection with string theory [2lJ , but is also of long-standing interest in the theory of non-linear waves 
[22j |. The Born-Infeld equation reads 



#tt(l + z x) ~ z xx{^ — z t ) — 2z x z t z xt 

Hoppe [23| gave a general solution of the form 

z'(t,<p) = Acos(/-ff)cos(/ + ff), z'(t,ip) = Xcos(f - g)sin(f + g) 
<P) = - sin (/ - .9) sin (/ + fiO, *(*) ¥>) = sin(/ - g) cos(/ + g), 



(75) 



(76) 
(77) 



where / = f((p + t/X) and g — g(ip — t/X) are any two (smooth) functions and if is a parameter. For (|76[) to be a 
graph, the tangent vector ([76]) should not be vertical, i.e. we must require that |/ + g\ < n/2. 
The curvature of this solution is 



f' + g' 

cos(/ - g) ' 



(78) 



so a singularity is expected whenever / — g = ir/2. Let us assume for simplicity that 5(C) = — /(— C)- I don't expect 
this to be a restriction on the class of possible singularities that can occur. Let f(Q have the local expansion 



/(C) = tt/4 + a(( - Co) - b(( - Co) 2 + 0(C ~ Co) 3 , 



(79) 
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so together with the symmetry requirement we find 

/ - g = tt/2 + 2a(t - Co) - 2% 2 + (t - Co) 2 ). (80) 

The factor A can be absorbed into a. 

From this expression it is clear that Co has to be identified with the singular time to and a > for the solution to be 
regular for t < to- To expand around the singular time, we put t' = to — t. Similarly, one must have b > (otherwise 
/ — g would be tt/2 at an earlier time), and the singularity occurs at (p = 0. Thus to leading order we have 

/ - g » tt/2 - 2at' - 2btp 2 , f + g^2ay, (81) 

from which we get 

x' = 2at' + 2b(p 2 , z' = 4a 2 t'tp + 4ab(p 3 . (82) 
Integrating this expression, using the integrability condition (1771) . gives 

x = t'lf + 2cip 3 /3, z = iV/2 + c<^ 4 , (83) 
where we have used a rescaling of the parameter ip. This is of course exactly the swallowtail ([3]) with 7 = 1. 
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